Robustness of spin-coupling distributions for perfect quantum state transfer 
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The transmission of quantum information between different parts of a quantum computer is of 
fundamental importance. Spin chains have been proposed as quantum channels for transferring 
information. Different configurations for the spin couplings were proposed in order to optimize the 
transfer. As imperfections in the creation of these specific spin-coupling distributions can never be 
completely avoided, it is important to find out which systems are optimally suited for information 
transfer by assessing their robustness against imperfections or disturbances. We analyze different 
spin coupling distributions of spin chain channels designed for perfect quantum state transfer. In 
particular, we study the transfer of an initial state from one end of the chain to the other end. We 
quantify the robustness of different coupling distributions against perturbations and we relate it to 
the properties of the energy eigenstates and eigenvalues. We find that the localization properties of 
the systems play an important role for robust quantum state transfer. 

PACS numbers: 03.67.Hk, 03.65.Yz, 75.10.Pq, 75.40.Gb 



I. INTRODUCTION 

Quantum information processing has been extensively 
studied during the past years [ ]. One of the main chal- 
lenges of actual physical implementations has been the 
manipulation of the quantum information with sufficient 
accuracy to prevent errors. In particular it is important 
to be able to transfer quantum information between dif- 
ferent elements of a quantum computer [2]. In this re- 
spect, spin chain systems have been proposed as quantum 
channels for the transmission of quantum states, where 
the spins act as the quantum bits [3-10]. Many systems 
of this kind have been explored in order to improve their 
performance for the state transmission. One of the goals 
is to find systems that allow for state transfer without 
any dynamical manipulations during the transfer proce- 
dure or with only minimal additional requirements. For 
example, local control only on the boundary spins in ei- 
ther an initialized chain [11, 12] or an unpolarized chain 
[13, 14] can cause a large enhancement of the transmis- 
sion fidelity from one end of the chain to the opposite 
end; even perfect state transfer (PST) could be achieved 
by engineering the entire set of spin-spin couplings in the 
chain [15-18]. 

Very few of these systems have been implemented ex- 
perimentally, for example using small numbers of spins 
in liquid state NMR [19-23] and slightly larger numbers 
of them in solid-state NMR [24, 25]. Spin defects in dia- 
mond seem to show a promising direction for near future 
implementations [26—28]. Important experimental chal- 
lenges are posed by the lack of individual addressibility 
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of the spins and, more importantly, by their vulnerabil- 
ity to decoherence [29]. Imperfections in the implemen- 
tation of spin-chain systems also cause decoherence and 
were predicted to produce localization of the quantum 
information [30-33] which was recently demonstrated ex- 
perimentally [34]. Consequently, a successful character- 
ization of these PST protocols should consider these er- 
rors in order to find the optimal one. Two PST proto- 
cols that require engineered spin-couplings [16, 17] have 
been analyzed in this respect considering static perturba- 
tions [18, 30, 35]. Other important points to consider are 
the timing errors on the time when the PST is achieved 
[18] and the speed of transfer of the different protocols 
[36]. But, considering that the number of possible sys- 
tems that could be used for PST [17, 18] is infinite, a 
performance comparison between them should be aimed 
at finding a system which is as robust against pertur- 
bations as possible. For that purpose it is important to 
find out which intrinsic properties of a system make it 
robust against perturbations. In this work, we tackle 
these questions analyzing different energy distributions 
that allow for PST and compare their robustness against 
static perturbations. We characterize the robustness of 
the systems by calculating their transmission fidelity. In 
order to find the relevant properties of the systems that 
make them robust, we analyze how the eigenstates and 
eigenenergies are perturbed. We find that the localiza- 
tion properties of a system are intimately connected to 
its robustness. 

The paper is organized as follows, in Sec. II we present 
the XX model describing the quantum spin chain and the 
necessary conditions for perfect state transfer. In Sec. Ill 
we analyze different energy eigenvalue configurations of 
the system and the corresponding spin-coupling distribu- 
tions. In Sec. IV A we analyze the fidelity of the trans- 
fer of the different configurations, and the influence of 
perturbations on the transmission is discussed in Section 
IV B. Subsequently, in Sec. V, we analyse how the indi- 
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vidual perturbed eigenstates and eigenvalues contribute 
to the dynamics of quantum information transport. Fi- 
nally, in Sec. VI we give the conclusions. 

II. PERFECT STATE TRANSFER CHANNELS 

We consider a chain of N spins 1/2 (qubits) with a 
modulated XX interaction between nearest neighbors. 
Taking into account an external magnetic field, the 
Hamiltonian is 

H = E 7W1 + - E |°f . 

i i 

where af are the Pauli operators of the ith spin, bi is the 
local external field and Jj is the exchange coupling. 

The aim is to transmit a quantum state |^o) stored on 
the first spin (i = 1) to the last spin of the chain (i = N), 
where |^o) — ot |0) + /3 |1) is a given superposition of a 
spin down and up respectively and the remaining spins 
of the chain are initialized in the spin down state. The 
Hamiltonian (1) preserves the total magnetization along 
the z-axis because [H, Sjaf] = 0, i.e., the number of 
excited spins is conserved. Because the initial state \^q) 
is a superposition of the eigenstate |0) = |00...0) and the 

state |1) = 1 1 1 0), the component |0) is conserved 

and the component |1) evolves within the one excitation 
subspace spanned by the basis states |i) = |0...01,0...0). 
The state of the system at a given evolution time t is 

N 

|*(t)> = e-^/ fi |* ) = «|0) + /*(*)!*>> ( 2 ) 

»=i 

where /j(t) = {i\e~' lHt / n \l). To measure the effectiveness 
of state transfer between sites 1 and N, we determine the 
fidelity F(t) = (^o\pN{t)\^o) averaged over all possible 
initial states |^o) distributed uniformly over the Bloch 
sphere, which is given by [3] 

I f N (t) I cos 7 , \f N (t)\ 2 1 , . 

F(t) = 3 + —^+2' (3) 

where 7 = arg |/jv(t)|. Because the phase 7 can be con- 
trolled by an external field once the state is transferred, 
we consider cos 7 = 1. PST is achieved when F = 1. 

For a spin chain possessing mirror symmetry with re- 
spect to the center, i.e., Jf = Jfj_ i and bi = b^+i-i, the 
necessary and sufficient condition for PST is 

Wfc+i - oj k = (2m fe + l)7r/ipsT, (4) 

where the set of eigenenergies {cok} is ordered, uj^ < 
uik+i- The condition (4) must be fulfilled for all pairs 
of successive energies, where the nik may be arbitrary 
integers. The shortest time ipsT for which (4) is fulfilled 
is the first time at which PST is achieved [17, 37]. Since 
(4) implies strictly periodic time evolution, PST occurs 
again and again, at all odd multiples of tpsT- 



III. ENERGY AND SPIN-COUPLING 
DISTRIBUTIONS 

Every set of integers mk in (4) leads to a unique energy 
spectrum enabling PST and hence, as we shall explain be- 
low, to a unique set of coupling constants Jj . Therefore, 
there are infinitely many spin chains allowing for PST. 
But, are all of them equally efficient for transferring in- 
formation? How is their PST capability affected by per- 
turbations through inaccuracies in the coupling constants 
or from coupling to external degrees of freedom? What 
properties are necessary to stabilize the system against 
such perturbations? 

We tackle these questions by studying the transmission 
robustness of different PST channels in the presence of 
static perturbations. We characterize these spin-channel 
systems by their energy eigenvalue distributions. A given 
spectrum that satisfies the condition (4) defines a unique 
Hamiltonian with positive symmetric couplings Jj , which 
can be obtained by solving an inverse eigenvalue problem 
[38]. For simplicity we choose Wfc = — u>N+i—ki Vfe, which 
imposes bi = 0, Vi [37]. 

In order to study a range of different eigenvalue dis- 
tributions systematically, we start from the case of an 
equidistant energy spectrum, nik = const, in (4), which 
was discussed in Ref. [16]. We vary that spectrum by 
distributing the energy values more densely either in the 
center or towards the boundaries of the energy spec- 
trum. The class of energy spectra which we discuss can 
be parametrized as follows: 

Uk{k/3,a) = -A/ k a )Sgn(k - k )[(kp -\k- k \) a - kg]. 

(5) 

We assume that is odd, k = l,...,N numbers the 
energy eigenvalues in ascending order, as before, and 
ko = jV 2 ~ 1 marks the center of the spectrum. The 
shape of the spectrum is controlled by an exponent a 
and a reference index kp which can assume two values; 
kp = kb = ko or kp = k c = 0. The overall width of the 
spectrum is controlled by Q ) . The equidistant energy 
spectrum (constant density of eigenvalues) is given by 
ojk(k c , 1). The density of eigenvalues in the center of the 
spectrum increases for both u)k(k c , n) and Wk(kb, — ) with 
integer n > 2. A larger density of eigenvalues close to 
the boundaries of the spectrum is obtained for Wfc(fcb,n) 
and u!k(k c , -). The shapes of the two spectra defined 
by these two possibilities for given n are different, as 
are those of ojfc(fc c ,n,) and ui^kb, —), respectively. For 
non-integer exponent a the energies (5) normally do not 
fulfill the commensurability condition (4) and have to be 
slightly readjusted to make PST possible. Figure 1(a) 
shows the energy eigenvalues for the equidistant spec- 
trum, Wfc(fc c , 1), along with the four possibilities just dis- 
cussed, for n = 2. The corresponding exchange couplings 
Ji (normalized by the maximum coupling strength J max ) 
are shown in Fig. 1(b). The coupling distribution deter- 
mines the transmission velocity as we shall discuss in Sec. 
IV. 
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Figure 1: (Color online) (a) Energy eigenvalue distributions 
LJk{kp, a). The symbols represent the energy values and the 
lines give the exact functional dependence of Eq. (5). (b) Ex- 
change couplings determined by solving the inverse eigenvalue 
problem for each of the spectra given in (a). 



IV. PERFECT STATE TRANSFER STABILITY 
OF ENERGY DISTRIBUTIONS 



A. Unperturbed transfer 



Fig. 2(a) and (b) show that the linear Wfc(/c c ,l) and 
quadratic u>k(k c ,2) distributions achieve perfect transfer 
without secondary maximum of the fidelity at some ear- 
lier time. For those two systems tpsr is thus equal to 
the ty[ of the spin chain. To make a quantitative analysis 
of the speed of the transfer, we consider as a reference 
the known value t\j for a homogeneous spin-chain with 
J L = J in Eq. (1). In that system, constructive interfer- 



ence at site N occurs at time t 1 ^ 



[42]. The transfer 



obtained at that instant is not perfect, but by switching 
couplings on and off to perform consecutive swap oper- 
ations, perfect transfer may be achieved at t s ^ ap ~ ^ 
[43]. In terms of the maximum coupling J m ax, the PST 
time for the linear distribution is t l p r g e p r = 4 J rJV , which 
is two times faster than the consecutive swaps assuming 
J = J m „v, t l %% e rP r = ht s ™ ap , but slower than the free 



j-linear 
■ l PST 

evolution, t l pg e p r = 



,sioap 

in a homogeneous chain. The 
other distributions are about 15 times slower than the 
linear case as listed in the caption of 2. The second im- 
portant time scale is given by the duration of the PST 
maximum of the fidelity, i.e. the time during which the 
fidelity is very close to unity. At can be interpreted as 
the time of residence of the perfectly transmitted state 
on the last site of the chain; it determines the timing 
precision required for perfect state read-out. While the 
quadratic distribution is much slower than the linear one 
in terms of transfer time, its advantage is a much longer 
window time. We will return to this point later in Sec. 
V. 



To compare the perfect state transfer performance of 
the spin-channels with the different energy eigenvalue dis- 
tributions of Fig. 1, we calculated their averaged fidelity 
(3). Figure 2 shows the fidelity of state transfer from one 
end of the chain to the other, as a function of time. The 
time scale is given in units of the first perfect state trans- 
fer time tpsT- At this point it is important to note that 
the dynamics of the system contain at least two other rel- 
evant time scales besides the time fpsT which we shall use 
as a unit of time. The first such time scale is the spin- 
channel clock time 2tM, i-e., the characteristic time of 
the information propagation within the chain, also called 
mesoscopic echo time [39]. For a chain supporting spin 
waves as elementary excitations, e.g. the uniformly cou- 
pled XX or Heisenberg chains, t^j is fixed by the max- 
imum group velocity of the spin waves [3, 40-42]. The 
group velocity of excitations with dispersion u>(k) (where 
k now temporarily denotes the wave number) is given by 
v g = ^7. Unfortunately this concept breaks down for 
the systems of interest here, since translational invari- 
ance is broken by the non-uniform couplings Ji and the 
wave number is no longer defined. Our numerical results 
show that tpsT can be larger than £m, see for example, 
Fig. 2(c) and (d). At tu the excitations created at t = 
at site i = 1 interfere constructively but not perfectly 
at site i = N. Perfect interference occurs only later, at 
t = tpgT, after the excitations have travelled back and 
forth between the ends of the chain many more times. 
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Figure 2: Averaged fidelity of the state transfer in a N = 31 
spin chain for the different energy distributions shown in Fig. 
1 as a function of time. The linear uik(k c ,l) and quadratic 
u)k(k c ,2) distributions achieve the perfect transmission with 
the first echo, while the other cases achieve it after several 
echoes. In panels (d) and (e), the black regions are due to 
fast oscillations because of the coupling strength oscillations 
shown in Fig. lb. The transfer times are given by ipg^ r = 
and tpsr ~ 7*PST r with 7 = 15.4, 17, 15, 14.5 for a 

quadratic, LJk(kb, 2)1 Wfc(fc c , 5) an d 0Jk(kb, 2) distribution re- 
spectively. 
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B. Perturbed transfer 

So far, we have discussed the performances of different 
spin-channels without any external perturbation. How- 
ever, since the perfect engineering of all spin couplings 
is highly improbable, the study of the performance of 
different spin-coupling distributions under perturbations 
by flawed spin couplings becomes relevant. To study the 
robustness of the spin chains against perturbations we in- 
troduce static random spin-coupling imperfections quan- 
tified by Si 



Ji 



Ji(l + Si), 



(6) 



where each 8i is an independent uniformly distributed 
random variable in the interval [—£,/, £j]. Sj is a positive 
real number that characterizes the maximum perturba- 
tion strength relative to Ji. 

We calculate numerically the fidelity time evolution 
F(t) = (F(t)) N averaged over A av different realizations 
of the random imperfection values Si. Figure 3 shows 
the averaged fidelity evolution for the different energy 
eigenvalue distributions for a common ej value. Only 
two cases are strongly distinguished by their robustness 
against the perturbation: the linear distribution, which 
was already studied by De Chiara et. al. [30], and the 
quadratic distribution. As the near-perfect echoes in Fig. 
3(a) and (b) show, disorder at the level of ej = iCP 2 
does not significantly affect PST in those two cases. In 
contrast, panels (c) and especially (d) and (e) of the same 
figure show a rather rapid decay of the fidelity (black line) 
down to a useless level. The colored lines in panels (c), 
(d) and (e) show the fidelities of the unperturbed systems 
for comparison. 
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Figure 3: (Color online) Averaged fidelity of the state transfer 
in a N — 31 spin chain with random perturbations of strength 
ej — 10~ 2 averaged over N av = 10 2 realizations for the dif- 
ferent energy distributions from Fig. 1 as a function of time. 
The colored lines in (c), (d) and (d) show F(t) for the systems 
without perturbation [sj = 0). 



Figure 4 shows in detail the comparison of the 9th echo 
between the linear and the quadratic distributions. The 
quadratic distribution is obviously more robust than the 
linear distribution, and also its At is larger. Also shown 
are the 9th echoes for n = 3 (colored line) and for n = 1 0, 
respectively. These data show that both the maximum 
fidelity and the length At of the time window for the 
state read-out increase with n for energy eigenvalue dis- 
tributions of type ujk{k Cl n). However, the increase from 
n = 3 to n = 1 is insignificant compared to the increase 
from n = 2 to n = 3. 
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Figure 4: (Color online) Averaged fidelity of the state transfer 
in a N = 31 spin chain with random perturbations of strength 
ej = 10~ 2 averaged over N av = 10 2 realizations for eigenvalue 
distributions ojfc(fc c ,n). Shown is the range of times around 
the 9th echo of the PST in the unperturbed chain, for n — 
1,2, 10 (black lines), and n — 3 (faint colored line very close 
to the n = 10 results). 



The increase of At with growing n can be explained by 
the changes of the exchange couplings Ji shown in Fig. 
1(b). When n changes from I to 2 the Ji decrease close 
to the boundaries and increase in the center of the chain. 
This trend continues even more strongly for larger values 
of n (data not shown) . The small spin couplings close to 
the boundaries of the chain prevent the spreading of the 
information once it is localized at one of the chain ends, 
thus leading to larger values of At. 

We now focus onn = 2 because the robustness proper- 
ties are similar for larger values of n, but it should be kept 
in mind that energy eigenvalue distributions ujk{k c ,n) 
with larger n are generally more robust. To determine 
the robustness of the different distributions, we calculate 
the averaged fidelity F(tpsT, £j) as a function of the per- 
turbation strength ej for the first PST time tpsT deter- 
mined from the unperturbed case. Figure 5 shows results 
for different energy distributions and for a wide range of 
perturbation strengths. The linear and quadratic dis- 
tributions turn out to be the most robust ones for all 
perturbation strengths of interest, yielding quite similar 
results for weak perturbations (ej < 0.2) where the fi- 
delity is larger than F = 0.9. For larger perturbation 
strengths, the quadratic distribution is most robust. 
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Figure 5: (Color online) Averaged fidelity at time tpsr as 
a function of the perturbation strength ej for the different 
energy distributions from Fig. 1 with N = 31 and N av = 10 2 . 



Recently, it has been shown that the relative decay 
of the mesoscopic echoes between a perturbed evolution 
and the corresponding unperturbed evolution could be 
used to determine and characterize the decoherence time 
of the spin-chain channel [44]. Similarly, to determine 
the decoherence time for each perturbation strength, 
we study the state transfer fidelity for different PST 
echoes as a function of their respective PST echo times 
tp ST = (2i — l)tpsT ; i-e., the times where the i-th PST 
echo arrives at site N for an unperturbed evolution. Fig- 
ure 6 shows the fidelity F (tp ST , e,j) as a function of tp ST , 
for different perturbation strengths ej. The left panel 
shows the fidelity for the quadratic distribution while the 
right panel illustrates the linear distribution. The deco- 
herence time, i.e., the decay time as a function of t P g T is 
longer for the quadratic distribution than for the linear 
one. Additionally the fidelity of the quadratic distribu- 
tion converges to an asymptotic value higher than that 
of the linear distribution due to the localization effects 
caused by the small couplings in the borders, as discussed 
above. Weak coupling between terminal qubits and the 
intervening spin chain were used as key elements also in 
other proposals for quantum information transfer by spin 
chains recently [11-14, 45]. 



V. ROBUSTNESS AND LOCALIZATION 

We have shown that certain systems are more robust 
against perturbations than others. In order to optimize 
the engineered spin coupling distributions it is decisive to 
understand which properties of the system are relevant 
for the robustness of the perfect state transfer. To this 
end we determine how each energy eigenstate contributes 
to the dynamics for each of the spin-channels. We ex- 
pand the states \i) (a single excitation at site i) in the 
eigenstate basis |i) = J^ fe dk,i\^k), where k numbers the 
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Figure 6: (Color online) Averaged fidelity at odd multiples of 
tpsr (symbols) for the linear and quadratic energy eigenvalue 
distributions, uik{k c ,n) (n = 1,2). Chain length is N = 31, 
averages were performed over N av = 10 2 realizations. Pertur- 
bation strengths are e = 0.01, 0.02, 0.1. 



energy eigenstates in ascending order, as usual. Figure 7 
shows the weights P^j = |afe.i| 2 , for the different energy 
spectra from Fig. 1. The mirror symmetries with respect 
to both the center of the chain and the center of the en- 
ergy spectrum are due to the spatial mirror symmetry 
of the couplings, and the symmetric tridiagonal nature 
(with zero diagonal) of the Hamiltonian matrix, respec- 
tively. Under perturbations of the couplings the spatial 
symmetry of the patterns of Fig. 7 is destroyed, while the 
energetic symmetry is not. We observe that the degree 
of localization of the energy eigenstates varies strongly 
between the different eigenvalue distributions. The most 
robust distributions seem to generate the most strongly 
localized energy eigenstates; in panels (a), (b) and (c) 
of Fig. 7 each energy eigenstate basically seems to be 
localized on two lattice sites. The quadratic distribu- 
tion [panel (a)] seems to have the most strongly localized 
eigenstates; in particular the eigenstates that belong to 
the center of the band are highly localized on the bound- 
aries. This is particularly clear from the upper panels 
in Fig. 7, showing the contributions Pk.i of the energy 
eigenstates \k) to the initial state |i = 1) with a sin- 
gle excitation localized at the boundary of the chain. In 
comparison, the other energy distributions show a larger 
spread in the contributions of the energy eigenstates to 
each site eigenstate \i). Nevertheless, we observe sim- 
ilarities of the distribution of Pk,i between the linear, 
quadratic and u>k(kb, h) distributions. It has been shown 
that the presence of localized states at the boundaries 
of the spin chain can improve the transmission of quan- 
tum states [12, 14, 46, 47]. These localized states arise 
when the coupling of the boundary sites is weaker than 
the coupling between inner sites or if external fields are 
applied at the boundary sites. Therefore, we study how 
the different energy levels are affected by perturbations 
for the different energy distributions. 
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Figure 7: Eigenvector probability Pk,i of the site (computational) states \i). Pkj = a\ iy where \i) = afc i j|\E'fc). The top part 
of each panel shows the probabilities P^j of the initial state |^o) = |1), and thus shows which energy eigenstates contribute to 
the state transfer. The panel labels refer to the different energy distributions given in Fig. 1, where a) u>k(k c ,2) quadratic, b) 
LUk(k c , 1) linear c) uj k (k b , |), d) ujk(k c , |), and e) ujk(k b , 2). 



We generated distributions of energy eigenvalues Wk for 
the different kinds of unperturbed energy spectra and for 
different perturbation strengths ej. For small Ej we ob- 
serve a symmetric distribution of the perturbed eigen- 
values u)k around their respective unperturbed values. 
The width of that distribution scales with the pertur- 
bation strength. For larger values of ej the distributions 
of the perturbed LOk become asymmetric with respect to 
the unperturbed energy level; the low-lying levels tend 
to be pulled down, while the high-lying levels are pushed 
up by the same amount. (The energy spectrum of the 
perturbed Hamiltonian matrix is still symmetric.) The 
value of Ej where the asymmetry sets in depends on the 
type of unperturbed energy spectrum and is largest for 
the quadratic case. To see more quantitatively what is 
going on in detail, we show in Fig. 8 the standard de- 
viations of the energy levels for the different kinds of 
unperturbed spectra. Each data point represents an av- 
erage over N av = 10 3 realizations of the random pertur- 
bations. The symmetry of the data with respect to the 
center uj = of the energy spectrum and the fact that 
the zero energy eigenvalue is not affected by the random- 
ness at all are due to the nature (symmetric, tridiago- 
nal, zero diagonal elements) of the Hamiltonian matrix. 
The key observation explaining the differences in state 
transfer robustness is made by combining figures 8 and 
the upper panels of Fig. 7. Those panels show that for 
the quadratic energy spectrum Wfc(fc c , 2) the initial state 
\i = 1) of the state transfer process is superposed from a 
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Figure 8: (Color online) Standard deviation of the energy lev- 
els |Aw ej | = \lo £ j —ujo\ due to the perturbation with strength 
ej for the different energy distributions of Fig. 1. For weak 
perturbations (small ej) the standard deviation turns out to 
be proportional to ejoj ma x, which we use as a unit here. The 
data shown are for ej < 0.1; N av — 10 3 realizations were used 
for the calculations. 



small number of energy eigenstates in the center of the 
energy spectrum. In all other types of energy spectrum 
the initial state shows a wider distribution in the energy 
quantum number k. At the same time, the sensitivity to 
perturbations (which is what is shown in Fig. 8) shows 
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a comparatively wide minimum with value zero in the 
center of the spectrum, whereas all other types of spec- 
tra show roughly constant nonzero values in the central 
region of the spectrum, with a single exceptional zero 
right in the center. The observed particular robustness 
of the quadratic energy spectrum can thus be ascribed 
to the fact that the initial state consists of a particularly 
small number of energy eigenstates coming from a part 
of the energy spectrum which is particularly insensitive 
to perturbations in the spin coupling constants. For the 
less robust distributions, the variance is roughly indepen- 
dent of the energy in the center of the energy band, while 
for the quadratic energy spectrum the variance decreases 
continuously towards the band center. For that distri- 
bution (and for the two other distributions shown in the 
upper panels of Fig. 7) all energy eigenstates are quite 
strongly localized. A glance at Fig. lb shows that the 
corresponding coupling patterns have the smallest cou- 
plings close to the ends of the chain, in the region where 
those energy eigenstates are localized which are most im- 
portant for the state transfer. Since we discuss a constant 
relative strength s of the disorder, the absolute changes 
of the couplings tend to be smallest near the ends of the 
chain, causing only small changes in the energy eigen- 
values. That explains the particular robustness of the 
quadratic distribution. 

Another important aspect characterizing the robust- 
ness of the transmission is the length of the window of 
time where high fidelity is obtained for the transmitted 
state. In this context, we are not only considering er- 
rors in the engineered spin couplings, but also the timing 
error of the measurement [18]. To that end we analyze 
the term \f N (t)\ 2 = |(N|e- ifft / R |l)| 2 from (3) at time 
tpsT + St. Taking into account the spatial symmetry, 
|/jv(*FST + St)\ 2 is given by 

I/aH 2 =|^(-l) fc+s P 5 ,i^,ie- l( " fc -" s)( * PST+5t) 

k . s 

«| P s,xPk,i(l - St{u k - u s )+ 

k . s 

St 2 

* 2 -^-( w fc - w s ) 2 - ...)| 

St 2 

k.s 

where P k .i = K^fcll)! 2 . Even without the Bloch-sphere 
average (3) which would be necessary for a comparison 
with Fig. 2, the above result shows why the quadratic 
distribution displays the longest window of time. The 
probabilities P k i (see Fig. 7) are sharply peaked in the 
center of the energy band and essentially zero otherwise. 
Furthermore, due to the quadratic nature of the energy 



spectrum, the relevant energy differences ui k — oj s are 
particularly small (see Fig. la) making the fidelity 
deviate from unity only at rather large St values. 

VI. CONCLUSION 

We have studied the robustness of spin chain systems 
designed for perfect quantum state transfer (PST) under 
static perturbations. We explored different PST systems 
by choosing different energy spectra distributions that 
satisfy the PST conditions. From the energy spectrum 
of a given chain, the spin-spin coupling constant pattern 
can be obtained by solving an inverse eigenvalue problem. 
The robustness of each system was studied by calculat- 
ing its transmission fidelity under static perturbations of 
the couplings. We found that robustness is characterized 
by two main features. One is the reduction of the trans- 
fer fidelity induced by the perturbed couplings and the 
other is the duration of the time window during which the 
transmitted state may be read out with high fidelity. The 
most robust systems are those with linear and quadratic 
energy eigenvalue distributions. These systems achieve 
PST at the time of the first fidelity maximum. That 
time may be called the spin-wave echo time, and the less 
robust systems reach PST only after several spin-wave 
echoes. By analyzing how the energy eigenstates and 
eigenvalues are affected by the perturbations, we found 
that the most robust distributions have strongly spatially 
localized eigenstates. Thus, because the initial state is 
localized in one end of the chain, only few eigenstates 
participate in the transfer. Because of the localization 
properties of the eigenstates the perturbations in the spin 
couplings close to the chain boundaries are the only sig- 
nificant source of errors. Since these couplings are rather 
weak for the most robust systems, a given relative pertur- 
bation strength only causes a small absolute perturbation 
in the couplings and thus in the energy eigenvalues, lead- 
ing to the observed robustness. The weak couplings close 
to the ends of the chain also lead to a longer residence 
time of the transmitted state at its target site at the chain 
boundary, causing a longer time window for read-out. 
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